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In turbulence modeling, the two-component limit is defined by (U3U3) = 0, or more 
generally by (uiUj)^j = 0 for some vector £, (here and throughout, Ui denotes ve- 
locity fluctuations). Although this limit certainly represents a rather extreme state of 
turbulence, it is approached in several practically important problems: for example, in 
near- wall turbulence, (U3U3) <C {U 2 U 2 ) ~ {uiUi). The two-component limit is also con- 
nected to the realizability problem: the condition that a predicted Reynolds stress tensor 
Rij can be represented as the correlation R^ = {uiufl for some real stochastic field u* 
is equivalent to the positivity condition Rij£i£j > 0 for arbitrary real indeterminates & 
(Schumann (1977)). If a stress field becomes unrealizable, so that Rij£*£j < 0 for some 
vector £* . it obviously must do so by passing through a two-component state in which 
Rij£i£j — 0- Formulating conditions on models that insure realizability a priori must 
therefore consider the behavior of models at or near the two-component limit (Shih and 
Lumley (1985)). Because of these considerations, the two-component limit has been given 
considerable attention in the modeling literature (Sjogren and Johansson (2000)). 

The purpose of this paper is to explore some wider implications of the two-component 
limit for both single-point turbulence models and spectral closure theories. Although the 
two-component limit arises most naturally in inhomogeneous problems like wall-bounded 
turbulence, the analysis will be restricted to homogeneous turbulence. But since homo- 
geneous turbulence is the crucial case for realizability (Girimaji (2004)), the conclusions 
will nevertheless be applicable to modeling. 

The essential point of our argument is that whereas the evolution of the stochastic 
velocity field is Markovian because it is governed by the Navier-Stokes equations, the 
exact stress evolution equation is not Markovian because it is unclosed. This property 
of moment evolution has been stressed by Kraichnan (1959). We will show that mod- 
eling stress evolution at the two-component limit with a closure that is Markovian in 
the stresses alone leads to basic inconsistencies in single-point modeling and, perhaps 
surprisingly, in spectral models as well. 

Turbulence transport modeling begins by assuming that the fluctuating velocity field 
can be described adequately by some finite number of basic statistics, for example, by 
the energy and dissipation rate in a two-equation model (Jones and Launder (1972)), 
or by the dissipation rate and the components of the stress tensor in a second moment 
closure (Hanjalic and Launder (1972)). Whilst such assumptions about the kinematics of 
turbulence are well satisfied in a very broad class of turbulent flows, turbulence modeling 
requires a much stronger assumption about dynamics, namely that these basic statistics 
satisfy closed equations of motion. However, no fundamental considerations justify this 
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assumption, which is instead supported after the fact by comparing model predictions 
with experimental measurements or numerical data. 

This paper proposes a critical assessment of this type of closure assumption in the 
important special case of second moment closure. We begin with the exact Reynolds 
stress transport equation in the form 

R-ij — Pij — — tij + Tij ( 0 . 1 ) 

The notation and terminology of (Girimaji (2004)) are used: the dot denotes the con- 
vective derivative with respect to the mean flow, Rij is the Reynolds stress, P \j is the 
production, is the pressure-strain correlation, dj the dissipation tensor, and T the 
transport. As noted earlier, in second moment closure, the basic statistics are the stress 
components and the dissipation rate; all other correlations will be expressed in terms of 
them and the mean velocity. The production is an explicit function of the stress itself 
and does not require closure. We will consider homogeneous turbulence, in which the 
transport term pj vanishes; thus, the pressure-strain correlation is the only unclosed 
term. 

Second moment closure assumes the functional dependence 

0 = f(R, VU, e) (0.2) 


where 0 denotes the pressure-strain correlation tensor, R the Reynolds stress tensor, 
and VU the mean velocity gradient: index-free notation has been introduced in order to 
emphasize that the left side is a tensor function of the arguments of f. This representation 
includes the possibility of a term of the form A(R, VU, e)l where I is the identity matrix 
and A is a function of the joint scalar invariants of R and VU and of e. 

We will show that Eq. (0.2) is inadequate in the two-component limit. The reason is 
that the pressure-strain correlation dynamics makes the evolution of the stresses non- 
Markovian. However, substitution of Eq. (0.2) in Eq. (0.1) results in a Markovian model 
for the stress evolution. Although a Markovian model of a non-Markovian physical phe- 
nomenon is never desirable, this simplification of the dynamics will prove to have par- 
ticularly dire implications in the two-component limit: in this limit, Eq. (0.2) yields a 
mathematically incomplete representation of the pressure-strain correlation leading to 
qualitatively inadmissable behavior, regardless of the functional form assumed for f. 

For this demonstration, we consider initially two-component turbulence for which 
R 33 (0) = R 13 (0) = R 23 (0) = 0. The exact equation for i? 33 is 




R33 2 <>«"r> ~ d'X3 / ~ \dXpdXp / 

and since all correlations on the right side vanish initially, also 

R 33 ( 0 ) =0 


(0.3) 


(0.4) 


Because every correlation in Eq. (0.3) contains a factor u 3 , it might appear that a per- 
manent two-component state in which w 3 = 0 is a solution. We will prove that instead, 
the Navier-Stokes equations drive the solution away from the two-component state. To 
rule out direct linear couplings that trivially generate nonzero w 3 , suppose that if a mean 
velocity gradient exists, it is nonzero only in the 1-2 plane. Then 


dU 3 _ dUj 

dxi 8x3 


(0.5) 


It will be convenient to say that the Reynolds stress and mean velocity gradient are 
two-component tensors. In this case, the production term on the right side of Eq. (0.3) 
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vanishes leaving 


R33 = —2 ( 113 


9p_\ / dus du3 \ 


8 x 3 / \ dx p dx p / 

The first correlation on the right side of Eq. (0.6) satisfies 

whereas the second correlation satisfies 

d_/^3 ^3 


dt \dx I 




Therefore, 


^ o)=2 (Z { *’ o) Z ix ’ o) 


The two-component limit certainly permits {(dp/dx 3 ) 2 ) |t=o 7 ^ 0. Therefore, 

^33(0) > 0 


(0.6) 

(0.7) 

(0.8) 

(0.9) 

(0.10) 


This result, with Eq. (0.4) implies that i? 33 (f) > 0 for sufficiently small positive t. Thus, 
the Navier-Stokes dynamics drive the solution away from the initial two-component state. 
Note that Eq. (0.10) states the ‘strong realizability’ condition (Shih and Lumley (1985), 
Girimaji (2004)). 

Now consider the predictions of second moment closure for this problem. Standard 
representation theorems for tensor functions (Pope (2000)) applied to Eq. (0.2) give 


R = A 0 I + AnQU (0.11) 

l<n<N 

where the C ra are an appropriate basis of tensors that can be formed from R and VU, and 
the A m with 0 <m<N are functions of scalar invariants of R and VU. In the case of no 
mean flow, the basis tensors C" consist of R and its powers. f All of the powers R" are two- 
component for n > 1 because R is, and Ao must vanish when t = 0 because i? 33 ( 0 ) = 0. 
It follows that R(0) must be two-component if there is no mean flow. If there is a nonzero 
mean gradient, the basis C" will also contain products of R and VU. However, since 
VU has been assumed to be two-component (Eq. (0.5)), all of these products are also 
two-component when t = 0. It again follows that R(0) is two-component. 

If R(0) is two-component, then 

i?13 (0) = i? 23 (0) = R33 (0) =0 (0.12) 

This conclusion is easily shown to be inconsistent with the equations of motion, which 
imply that 

i?is(0) = <«i(0)u 3 (0)) + ( Ul (0)u 3 (0)) (0.13) 

Although the first correlation vanishes, the second need not vanish in general. Thus, 
the assumed second moment closure form Eq. (0.2) is inconsistent with Navier-Stokes 
dynamics. This inconsistency leads to an even more serious problem. 

Consider the evolution of the ‘non- two-component’ stresses R13, R23, and R33 under 
Eq. (0.11). Because VU is two-component, some of the C” might happen to be two- 
component for arbitrary R. Such terms will not contribute to the evolution of any of i?i 3 , 


f In fact, the Cayley-Hamilton theorem implies that the explicit representation depends only 
on powers up to the second, but this fact is not essential to the argument. 
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i? 23 , or R 33 . But in general, the C" will not be identically two-component. Eq. (0.12) 
requires that each term A n C n be two-component when t = 0 for arbitrary values of 
i?n(0), i?i 2 ( 0 ), and i?i 3 ( 0 ). Consider the Taylor series expansion of A n C ^ for i = 1,2,3 
in powers of the stress components. Each term in the expansion must contain either i?i 3 , 
i? 23 , or i ?33 to at least the first power, otherwise A n C could only be two-component 
when t = 0 for special choices of i?n(0), i?i 2 ( 0 ), and i? 22 ( 0 ). We conclude that 

A n Cl = a{ n 3 R 13 + ajgiJaa + a^Rss for i = 1, 2, 3 (0.14) 

where the a m are functions of the components of R and VU. 

Substituting the results Eq. (0.14) in Eq. (0.11), it is seen that under the closure 
hypothesis Eq. (0.2), compliance with the initial conditions dictated by the equations of 
motion requires that the evolution equations for R\ 3 , R 23 , and R 33 have the special form 

Rl3 = ^ 11^13 + A 12 R 23 + Ai 3 R 33 
R23 = A21-R13 + A22R23 + A23R33 

R 33 = A 31 -R 13 + A 32 R 23 + A 33 R 33 (0.15) 

where the quantities Aij are functions of the components of R and VU and of e. Re- 
gardless of the analytical expressions for the Aij. Ri 3 (t ) = 0 satisfies Eq. (0.15) with the 
initial conditions i?* 3 ( 0 ) = 0; thus, the closure assumption of Eq. (0.2) implies that the 
solution to the model equations never leaves the two-component state. Thus, Eq. (0.2) is 
inconsistent with the equations of motion, which imply that the two-component state 
cannot persist. The problems of modeling near the two-component limit have been ana- 
lyzed by Craft and Launder (2002), who investigate the kinematic constraints imposed 
by the vanishing of certain components of the pressure-strain tensor at this limit. The 
present work complements this investigation by considering the dynamics of evolution 
away from a two-component state. 

The essential point can be illustrated in a simple model problem. Suppose that the 
stochastic variable x satisfies 

x = a (0.16) 

where a is a known stochastic process. Then 

j t {x 2 ) = 2 (ax) (0.17) 

Suppose that x(0) = 0. As in the stress equation, x = 0 is not a solution of the complete 
hierarchy of moment equations that begins with Eq. (0.17), because 

^(x 2 )| t=0 = 2 (ax)l t=0 + 2 (a 2 )l t=0 = 2 (a % =0 (0.18) 

need not vanish. But if we introduce a closure hypothesis that 

(ax) = F((x 2 )) (0.19) 

so that 

±(x 2 ) = F((x 2 )) (0.20) 

with F( 0) =0 imposed to respect the initial conditions, a differentiable F forces the 
conclusion that any solution for (x 2 ) that vanishes initially must in fact vanish identically. 
The problem is that since under Eq. (0.16), x ~ t for short times if x(0) = 0, (x 2 ) ~ t 2 
for short times, but (ax) ~ t. Since t is not a smooth function of t 2 near t = 0, the 



5 


Second moment closure near the two-component limit 

closure assumption in Eq. (0.19) cannot be correct. Arguing that, in view of the short- 
time behavior, we should set {ax) ~ \J {x 2 ) , the result is a singular differential equation 
(i d/dt){x 2 ) ~ \J (x 2 ), for which the solution satisfying (x 2 )(0) = 0 is not unique. 

The analysis for the stress tensor is essentially the same, although the presence of tensor 
components necessarily generates the more elaborate argument leading to Eq. (0.15). The 
common features are the Markovian evolution of x, Eq. (0.16), but the non-Markovian 
evolution of the moment (x 2 ) due to its unclosed evolution equation Eq. (0.17). The 
closure hypothesis Eq. (0.19), like Eq. (0.2), is inconsistent with the short-time behavior 
of {ax) and (x 2 ). 

This argument is restricted to the short-time properties of initially two-component 
turbulence. The demonstration that R33 grows is consistent with Kolmogorov’s ideas, 
according to which turbulence should become more isotropic as it evolves. We would like 
to supplement these short-time results with long-time results available from integrating 
the equations of rapid distortion theory (RDT) (Hunt and Carruthers (1990)). In the 
standard Fourier representation, these equations have the form 

Qij{k,t) =n ij (k,t) +§ij(k,t) (0.21) 

where 

Qij{k, t) = J drexp[ik-r]{ui{x + r,t)uj(x,t)) (0.22) 

is the Fourier transform of the two-point correlation tensor, 

njj(k,t) = Q ip (k,t)dU p /dxj + Qj p (k,t)dU p /dxi (0.23) 

is the production tensor, and 

$ij(k,f) = i{p{-k, t)[kjUi{k, t) +kiUj{k,t)]) (0.24) 

is the pressure-strain tensor in which the pressure p is evaluated in the linear approxi- 
mation as 

p{k,t) = k~ 2 k m u n {k,t)dU n /dx m (0.25) 

The Reynolds stress is expressed in terms of the correlation by 

Rij(t) = J dk Qij(k,t) (0.26) 

The important fact is that the RDT equations are closed evolution equations for 
Qij{k, t) because $jj(k,f) can be expressed explicitly in terms of Qij(k, t) through Eqs. 
(0.24) and (0.25). Of course, the integrated quantities Rij(t) and 

$ij(t) = J dk 4>jj(k,f) (0.27) 

do not satisfy closed equations; the ‘rapid pressure-strain’ model formulated as Eq. (0.2) 
attempts to express in terms of Rij, the mean gradients, and the dissipation rate. 
RDT is therefore a good vantage point from which to evaluate this hypothesis. 

Figure 1 shows the evolution of Rn{t),R 22 {t),R 33 {t), and R\ 2 {t) according to RDT, 
starting from a two-component <2ij(k,0) under mean shear dU\/dx 2 - The main con- 
clusion is that R33 (t) does indeed become nonzero; at large times, it is actually larger 
than R 22 • Figure 2 shows the evolution of $n(f),$22(f),$33(f), and $12 {t). This figure 
pinpoints the reason that Eq. (0.2) fails. Obviously $33 (t) ~ t at short times, whereas 
R 33 (t) ~ t 2 . These results can be compared to the model problem described earlier. As 
in that problem, second moment closure tries to express a quantity which grows linearly 
in time as a smooth function of a quantity which grows quadratically. 
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RDT shows that the coupling of the stress components to the rapid pressure-strain 
tensor sustains long-time growth of i? 33 (f). We would like to reconsider homogeneous 
anisotropic turbulence without mean shear from the viewpoint of spectral closure theory 
in order to understand the nonlinear mechanisms which lead to long-time growth of 
R 33 (t). 

In a classic paper, Herring (1974) studied the relaxation of axisymmetric turbulence 
using the direct interaction approximation (DIA) of Kraichnan (1959). We quote the 
relevant equations as 

<^> A (k; 1. 1') = —vk 2 (f> x { k; t, t') -1- ^ f dpdq <5(k — p — q) x 

n,v •* 

|A AjUJ/ ( k ,p,q) J ds g A (k;t,s)^(p;t,s)^(q;t,s) 

+ B x ^(k,p,q) ds g"(p-,t, s)<j> x (k;t, s)<f> v (q;t, s)^ (0.28) 

In this equation, A, yu, € {1,2}; they are not tensor indices, so any summations involv- 
ing them are indicated explicitly. ^> 1 (k;t, t') and (j) 2 (k]t, t') are scalar descriptors of an 
axisymmetric two-time correlation tensor in the elementary but elaborate kinematics of 
(Herring (1974)); isotropy is the special case q i 1 = cf> 2 with dependence on k = |k| only, 
and the two-component state corresponds to 4> 2 = 0. The quantities A Am „ and Bx^v 
are geometric factors. It is crucial to our argument that A AjU! , > 0; the H AjU! , are typi- 
cally negative, and represent damping effects. Finally, g x (k;t, s) are additional quantities 
linked to the temporal decorrelation of turbulence; we will only use the initial conditions 
3*(k; 0, 0) = 1. 

Herring (1974) shows by explicit calculation that initially two-component axisymmetric 
turbulence becomes fully three-component under the DIA evolution equations (see Figure 
9 of (Herring (1974)). The reason that DIA predicts the growth of (j} 2 (k;t,t') starting 
from ^> 2 (k;0,0) = 0 is that the A 2 n term in Eq. (0.28) couples ^ 2 (k ;f,f') to ^(k; t, t'). 
The short-time form of the evolution equation for cj) 2 (k;t, t') is 

4> 2 (k;t,t') = t' J dpdq d(k - p - q)A 2 n (k, p, q)^ M (p; 0, 0)<^(q; 0, 0) (0.29) 

Thus, (j} 2 (k;t,t') ~ tt 1 and the single-time correlation grows as <^> 2 (k,f) = (j) 2 (k,t,t) ~ t 2 ■ 
This short-time quadratic growth is consistent with earlier conclusions based on short- 
time properties of the Navier-Stokes equations. Note also that Eq. (0.29) shows that 
^ 2 (k, 0) > 0. 

We would like to highlight the ingredients in DIA which make departure from an initial 
two-component state possible: one way is to ask to what extent DIA can be simplified 
and yet remain consistent with this property. Suppose that we introduce the ‘almost- 
Markovian’ (Kraichnan (1971)) modification 


d(j) x 

dt 


(k, t) = —nk 2 <f> x (k, t ) + ^ f dpdq <5(k — p — q) 


{A AM „(k, p, q)0 A ^(k, p, q; t)^ M (p, t)4> v { q, t) 

- B x ^( k, p, q)0 MAy (k, p, q; t)(p x ( k, q, f)} 


(0.30) 


where <j> x (k,t) = (fi x {k,t,t) and the time-scales 0 A#tJ /(k, t) satisfy evolution equations 
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with the structure (Kraichnan (1971)) 

0 A ^(k,p,q;t) = 1 - [j? A (k,f) +VpiP,t) + 7?„(q, t)]0A^(k, p, q; t) (0.31) 

with suitable damping factors r]. The combination of Eqs. (0.30) and (0.31), which de- 
pends on single-time quantities alone, is a multiple state variable Markovian model like 
the test-field model (Kraichnan (1971)); Kraichnan calls such models ‘almost Markovian’ 
to distinguish them from models like EDQNM (Orszag (1973)) which are Markovian in 
the correlations alone. 

The short-time properties of this simpler model agree with those of Eq. (0.28): if, in a 
two-component initial state, we may assume that 0 Amj/ vanishes unless A = n = v = 1, 
then we will have 02ii ~ t for short times, and consequently 4> 2 ~ t 2 . But it is critical to 
this argument that ©211 evolve independently of 0m; the plausible attempt to simplify 
the model by using a single relaxation time 0 would lose these short-time properties. 
Despite its complexity, such a model would also predict the persistence of an initial 
two-component state. 

Another plausible but unsatisfactory simplification is to construct a theory that is 
Markovian in the correlations alone, for example, by replacing Eq. (0.31) by its long-time 
limit 0 AjUI ,(k,p,q;i) = l/[?j A (k, t) +rj^(p,t) +r] v (q, t)]. Regardless of how the damping 
factors rj are chosen, this theory cannot be made consistent with the short-time growth 
cf 2 ~ t 2 : given that the Navier Stokes equations will require ^> 2 (k, 0) = 0, this theory 
theory will force 4> 2 to vanish for all times. We conclude that (fully) Markovian spectral 
closure for axisymmetric turbulence is inconsistent with the equations of motion for two- 
component initial conditions. The minimal spectral closure model that is consistent with 
the properties of two-component turbulence is no more than ‘almost-Markovian.’ Time- 
scale evolution cannot be avoided, moreover, a complete matrix of time-scales must be 
evolved; a single scalar time-scale is insufficient. It should be noted that RDT can escape 
the two-component state using the correlation tensor alone because coupling to the mean 
field by the pressure-strain correlation makes the stress evolution non-Markovian. 

Although we have demonstrated that the closure hypothesis of Eq. (0.2) is incorrect 
only in a very special circumstance, this demonstration underscores that such hypotheses 
are uncontrolled approximations of uncertain validity. Perhaps more attention should be 
given to finding conditions under which these hypotheses can actually be justified. 

In conclusion, we would like to summarize some implications of this analysis for the re- 
alizability problem for single-point closures, which first motivated this investigation. We 
have found that at the two-component limit, the strong realizability condition i? 33 (0) > 0 
or its spectral analog: (k, 0) > 0 and ^ 2 (k, 0) > 0, is satisfied, whether as a conse- 

quence of short time analysis of single-point correlations, RDT, or DIA. To this extent, 
our results support the strong realizability constraint as a realistic statement of dy- 
namics near the two-component limit. However, we also find that the standard second 
moment closure assumption Eq. (0.2) is not consistent with this constraint, since it forces 
R 33 (0) = 0. It is significant that Markovianization of DIA also leads to the unsatisfactory 
conclusion that ^ 2 (k,0) = 0. 

The strong realizability constraint is therefore physically correct, but cannot be con- 
sistently enforced with the usual Markovian formulation of second moment closure. We 
believe that Markovianization of the time dependence underlies an important unresolved 
issue in rapid pressure-strain correlation modeling: fundamental constraints on these 
models are (i) that the closure be linear in the stresses, and (ii) that the closure satisfy 
the strong realizability condition. To date, no model satisfies both conditions, and it has 
been claimed that strong realizability requires nonlinear models. At the level of Marko- 
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vian second moment closure, the weak realizability condition (Pope (2000)), which makes 
the two-component limit inaccessible, therefore seems unavoidable. 

We would like to briefly outline some alternate approaches to the realizability problem 
suggested by the present analysis. One way to make a single-point model consistent with 
the two-component limit is to extend it, for example to include the components of the 
pressure-strain correlation as independent dynamic variables; for shear turbulence, the 
structure tensors (Kassinos et al. (2000)) are also reasonable candidates. This type of ex- 
tension of the model can be compared formally to Kraichnan’s almost-Markovianization: 
the introduction of a new variable abandons the Markovian character of Eq. (0.2) and 
allows R to be a functional of the time history of the arguments of f . 

Another possibility is reformulate the problem so that stress is not the basic variable. 
For example, (Girimaji (2004)) suggests that the realizability of the Reynolds stress can 
be guaranteed by ensuring the realizability of the closure models for each individual 
term in the exact stress transport equation. This approach shifts the enforcement of 
realizability away from the troublesome two-component limit. A related approach is to 
introduce a symmetric tensor S with the property that R = S 2 . If the theory is expressed 
near the two-component limit in terms of S rather than R, then a Markovian model for S 
will allow S ~ t near the two-component state, yet R ~ t 2 as required. The introduction 
of S effectively regularizes the problem near the two-component limit. It is critical that 
the theory is expressed in terms of S and not R: we are not suggesting the introduction of 
a matrix square root, which would lead to interminable problems of sign and the strong 
possibility of non-uniqueness in the time evolution. 
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